Mon. Not. R. Astron. Soc. 000, [T]-?? () 



Printed 1 December 2008 



(MN I^TeX stylo file v2.2) 



Molecular gas and star formation in the red-sequence 
counter-rotating disc galaxy NGC 4550 

00 . Alison F. Crocker/ Hyunjin Jeong,^ Shinya Komugi,^'^ Francoise Combes,^ 
^ ; Martin Bureau,^ Lisa M. Young, ^ Sukyoung Yi^ 

r ^ Sub- Department of Astrophysics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford 0X1 3RH 

' ^Department of Astronomy, Yonsei University, Seoul 120-749, Korea 

^ ' '^Institute of Astronomy, University of Tokyo, 2-21-1 Osawa, Mitaka-shi, Tokyo 181-8588, Japan 

O , National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka-shi, Tokyo 181-8588, Japan 

1 ^^ ^ ■ ^ Observatoire de Paris, LERMA, 61 Av. de I'Observatoire, 75014, Paris, France 

' ^Physics Department, New Mexico Institute of Mining and Technology, Socorro, NM 87801, U.S.A. 

o 

C^; ABSTRACT 

Q , We present observations of the CO(l-O) emission in the central 750 pc (10") of the 

$-H ■ counter-rotating disc galaxy NGC 4550, obtained at the Institut de Radioastronomie 

I Millimetrique (IRAM) Plateau de Bure Interferometer. Very little molecular gas is 

. detected, only 1 x 10^ Mq, and its distribution is lopsided, with twice as much molec- 
ular gas observed at positive relative velocities than at negative relative velocities. The 

I velocity gradient in the CO(l-O) emission shows that the molecular gas rotates like the 

^ i thicker of the two stellar discs, which is an unexpected alignment of rotations if the 

' thinner disc was formed by a major gas accretion event. However, a simulation shows 

, that the gas rotating like the thicker disc naturally results from the coplanar merger 

T-H ■ of two counter-rotating disc galaxies, demonstrating the feasibility of this scenario for 

' the formation of NGC 4550. We investigate various star formation tracers to determine 

, whether the molecular gas in NGC 4550 is currently forming stars. UV imaging data 

• and optical absorption linestrengths both suggest a recent star formation episode; the 

OO ! best-fitting two population model to the UV-optical colours yields a mass of young 

. stars of 5.9 x 10^ Mq with an age of 280 Myr. The best information on the current 

^ I star formation rate is a far infrared-based upper limit of only 0.02 Mq yr~^. We are 

• i-H . thus witnessing NGC 4550 either in a dip within a bursty star formation period or 

' during a more continuous low-level star formation episode. 
H ■ 

' Key words: galaxies: individual: NGC4550 - galaxies: elliptical and lenticular, cD 

- galaxies: ISM - ultraviolet: galaxies - galaxies: kinematics and dynamics - galaxies: 
stellar content 



1 INTRODUCTION 



Deciphering star formation histories of early-type galaxies is 
important for understanding their formation and evolution. 
While early-type galaxies were traditionally considered to be 
simple single stellar populations, it is now clear that their 
star formation histories must be more complica ted. Both 
absorption line strengths (e.g. iTrager et al.|[200ol) and UV- 
optical colours l|Yi et al.ll2005l : Kavirai et al. I 20071 ') indicate 
that a significant portion of early-type galaxies have recently 
formed stars (~ 30% from the UV-optical colours). Reassur- 
ingly, molecular gas, the raw mat erial for star formation, i s 
also found in early- type galaxies l|Henkel fc Wiklin"dlll997h . 
Two recent surveys give detection rates of 28% for E/SOs 
in the SAURON sample (|Combes. Young fc Bureaul 120071 ') 



and 78% for SOs in a volume-limited sample (ISage fc WelchI 
[2006). 

NGC 4550 is an unusual galaxy (see Table 1), containing 
two coplanar counter-rotating st ellar exponential discs with 
nearl y identical scale lengths (iRubin. Graham fc KennevI 
1 19921 : iRix et all 1 19921 '). However, integral field unit data 
shows that one disc is thicker and ha s a higher velocity dis- 
persi on than the other (|Kennev et al. 2000: Ems ellem et al.l 
l2004h . Schwarzschild modelling by [Cappellari et al.l (|2007l ) 
confirms the difference in scale heights and additionally 
shows that the two discs have equal mass within the 
SAURON field-of-view. The presence of these two discs with 
different scale height s may cause the failure of some bulge- 
disc decompositions. iRix et abl (| l992l) find a bulge to disc 
ratio of only 0.19 whil e [Scorza et al.l (Il998l) report a ratio 
of 5. The Sersic fit of iFerrarese et al.l l|2006h gives a best- 
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fit Sersic shape parameter of 1.7, demonstrating tfiat tlie 
surface briglitness profile is not dominated by a de Vau- 
couleurs bulge. With an effective |B — V| colour of 0.890 mag, 
NGC 4550 is very likely a red sequence galaxy, which would 
make it one of a very small population of disc-dominated 
red-sequence galaxies. 

Two main scenarios exist for explaining the cou nter- 
rotating stellar discs in NGC 4550. iRubin et al] (|l992l ') first 
suggested that accretion of counter-rotating gas and sub- 
sequent star formation could create a second disc rotating 
as observ ed. This gas accr e tion s cenario has been investi- 
gated by iThakar fc RvdenI l|l998h . who find it difficuH to 
create extended, exponential gas discs in their simulations. 
However, recent observations of NGC 5719 show an ex- 
tended counter-rotating a nd star-forming gas disc clearly ac- 
creted from external gas (|Vergani et al.l2007l ). demonstrat- 



Table 1. Basic properties of NGC 4550. The left and middle 
columns list the different quantities and their values; the right 
column lists the corresponding references. 



ing the feasibility of this origin. The second scenario invokes 
a coplanar major merger of two counter-rotating disc galax- 
ies. While requiring a precise alignment of the two discs 
to avoid excessive heating, this merger scenario has been 
shown to prod uce resultant galaxies resembling NGC 4550 
(jPuerari fc Pfe nnigcr 2001). 

The ionised gas in NGC 4550 co-rota tes with the thicker 
disc iRubin et al.ll 19921 : [Sarzi et al.ll2006l '). not with the thin- 
ner, colder disc. An irregular dust distribution is also ob- 
served within the central 20" in diameter. The dust is dis- 
tributed in clum py arcs and is stronger on the northern 
half of the galaxy (|Wiklind fc HenkefcOOll ) . HI observations 
have not detected neutral gas in NG C 4550, the strictest 
upper limit being 7 x lO'^ Mq from iDuprie fc Schneider] 
ll 19961). Molecular g as wa s first discovered in NGC 4550 by 
IWikhnd fc Henkell l|200ll ) using the IRAM 30m telescope. 
They reported a small mass of molecular gas, 1.3 x 10^ Mq, 
noting that the distribution of the molecular gas is likely 
to resemble that of the dust, due to a strong asymmetry 
in the observed CO emission line toward positive relative 
velocities. 

Here we present interferometric observations of the CO 
emission in NGC 4550, deriving the extent, kinematics and 
total mass of the molecular gas. To investigate what the 
kinematics of the gas reveal about the formation history of 
this unusual galaxy, we run a major merger simulation (im- 
portantly including gas) and describe the features of its evo- 
lution, comparing the results to observations. We also anal- 
yse various star formation indicators to determine whether 
the molecular gas in NGC 4550 alerts us to another early- 
type galaxy with ongoing star formation, or whether the 
molecular gas instead appears to be stable. 



Quantity 


Value 


Rcf. 


R.A. (J2000.0) 


12'' 35'" 30.6" 


1 


Dec. {J2000.0) 


-1-12° 13' 15" 


1 


Heliocentric velocity 


435 km s-i 


2 


Distance 


15.5 Mpc 


3 


Scale 


1" = 75 pc 


3 


Type 


SBO 


4 


Corrected apparent B mag 


12.31 


5 


Corrected absolute B mag 


-18.64 


5 


(B-V), 


0.89 


5 


Lb 


4.3 X IQS Lq 


6 


LpiR 


5.8 X 10'' L0 


6 


Lfir/Lb 


1.4 X 10-2 


6 


Lfir/Mhj 


8.0 L0/M0 


6 



References: ( 1) NED; (2) Derived from SAURON s tellar kine- 
matic data llEmsellem et al.l 120041') (3) iMei et al.l l|20o3); (4) 
Ide Vaucouleurs et al .1 lll991^ : (5) H yperLEDA; (6 ) Derived quan- 
tity using data from NED, IRAS l lMoshir et al.lll990l) and this 
paper. 



Table 2. CO observation calibrators 
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3C273 


Phase 


3C273, 1156+295 


Flux 


3C273, 1156-1-295, 0528+134, 3C84 
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2 PLATEAU DE BURE OBSERVATIONS 
2.1 Observations 

NGC 4550 was observed in the CO (1-0) line with the 
IRAM Plateau de Bure Interferometer (PdBI) on August 
10, 2007. The observations were taken in the D configu- 
ration with 5 antennas and used the new generation dual- 
polarization receivers at 3 mm. The spectral correlators were 
centred at 115.12 GHz, the transition frequency of CO(l-O) 
roughly corrected for the heliocentric velocity of NGC 4550. 
The correlator configuration used four 320 MHz-wide (834 



Figure 1. Channel maps of NGC 4550. The channels are 20 km 
s" ^ wide and contours are plotted at -3 (dashed), 3, 6, and 9 
a with (T = 2.77 mjy beam" ^. The number in the top left cor- 
ner of each frame is the central velocity of that frame, relative 
to the observed central velocity of NGC 4550 (14ys = 435 km 
s~^, determined using optical absorption lines). The synthesized 
beam is shown in the bottom left corner of each frame. The cross 
represents the centre of the galaxy as given by 2MASS. 
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Figure 2. Left: CO(l-O) integrated intensity map of NGC 4550. Right: CO(l-O) mean velocity map relative to the systemic velocity of 
435 km s~^ (see Table 1). The synthesized beam is shown in the bottom-left corner of each map. The cross represents the centre of the 
galaxy as given by 2MASS. 



km s~^) units with a frequency resolution of 2.5 MHz (6.6 
km s~^), covering a total usable bandwidth of 950 MHz 
(2475 km s~^). The correlator was regularly calibrated by a 
noise source inserted in the IF system. 

We obtained visibilities with series of thirty 45 s integra- 
tions on source, followed by three 45 s phase and amplitude 
calibrations on each of the two phase calibrators (see Table 
2). To flux calibrate, we set the flux density of 3C273 to the 
expected value (15 Jy), then checked that the fluxes found 
for the other flux calibrators were reasonable. The uncer- 
tainty in our flux calibration is ~ 20%. 

The data were reduced with the Grenoble Image and 
Line Data An alysis System (gildas) soft ware packages CLIC 
and MAPPING (|Guilloteau fc Luca^bOOd ) . We cahbrated the 
data using the standard pipeline in CLiC. After calibra- 
tion, we used mapping to create data cubes centred at 
115.104 GHz (corrected more precisely for a systemic ve- 
locity of 435 km s^^) with velocity planes separated by 20 
km s~^. Natural weighting was used. The primary beam size 
is 44" for CO (1-0) and we chose the spatial dimensions of 
the datacube to be about twice the diameter of the primary 
beam, 90" x 90". The synthesized beam is 4'.'7 x 4'.'0, so 
we choose spatial pixels of 1" x 1". Continuum subtraction 
was unnecessary (see below). The dirty beam has small side 
lobes that necessitated cleaning the dat acube. The clea ning 
was done using the Hogbom method l|Hogboml Il974h : we 
stopped cleaning in each velocity plane after the brightest 
residual pixel had a value lower than the rms noise of the 
uncleaned datacube. 

To constrain the continuum emission, we selected fre- 
quencies at least 40 km s~^ away from the lowest and highest 
velocity channels with any line emission in the cleaned dat- 
acube (-80 and 80 km s~^, respectively). The very edges of 
the bandwidth were also avoided, resulting in a 683 MHz 
wide continuum window. We mapped this data with the 
same spatial parameters as the line data. No emission could 
be seen in the map with an rms noise of 0.54 mjy beam~^, 
giving a 3a upper limit of 1.62 mJy assuming a point source. 



2.2 Results 

The channel maps (Fig. [TJ show CO (1-0) emission around 
the galaxy centre. The emission is just over the 3a level, 
with stronger emission in the 20, 40 and 60 km s~^ channels. 
The channel velocity values are all given with respect to sys- 
temic velocity. In the channel maps, a mild velocity gradient 
is observed, with more of the emission at positive relative 
velocities to the north of the galaxy centre and more of the 
emission at negative relative velocities to the south. The 60 
km s~^ channel is an exception with a peak slightly south of 
the galaxy centre, although the emission clearly also extends 
to the north. To pick up more extended emission, we created 
an integrated intensity map and a mean velocity field map 
(Fig. ^ using the smoothed-mask method. To make these 
maps, we first created a smoothed cube by smoothing with 
a 2D circular gaussian spatially (FWHM of 4" or about one 
synthesized beam width). Then we hanning smoothed by 3 
channels in velocity. The moments were computed only us- 
ing pixels in the original cube that corresponded to pixels 
above 3a in the smoothed cube. 

The integrated intensity map shows that the emission 
is limited to the central 750 pc (10") in diameter and is 
slightly stronger in the north. The mean velocity map shows 
a north-south gradient, as expected from the shift observed 
in the channel maps. Assuming the velocity gradient indi- 
cates rotation around the centre of the galaxy, the molecu- 
lar gas rotates like the observed ionized gas and the thicker 
disc. Spatially integrating the entire datacube over the re- 
gion with observed emission results in the spectrum shown 
in Fig. O This spectrum shows that the emission is present 
in all channels between -80 and 80 km s~^, but is strictly 
limited to these channels. 

Wh ile the lopsidodnoss of the CO emission was discov- 
ered by IWiklin d & Hcnkcl (2001), the de gree of asymme- 
try was exaggerated due to their use of an incorrect sys- 
temic velocity. The systemic velocity listed in NASA Extra- 
galactic Database ( NED ) for NGC 4550 and that used by 
IWikfind fc Henkell (|200ll i is 381 ± 9 km s'V This value is 



4 A. F. Crocker et al. 



based on the HI measurement of iPetersonI (|l979l ). More re- 
cent and more sensitive HI observations have not detected 
HI (DuPrie & Schnei der, 1996; M organti et al., in prepara- 
tion), making the iPetersonl l|l979i) value very questionable. 
In addition, recent absorption line stud ies seem to be con - 
verging to a value of around 435 km s ~^ l|Rubin et al.l (|l997l ) 
Falco et all (|l999l ) Hst 458 ± 41 km s"\ 
and the veloc- 



find 434 km s ., . . _ 

IWegner et al] (|2003l ) find 437 ± 15 km s 



ity in the central pixel of the SAURON velocity map is 
435 km s~^). We have adopted this value for our work. Then, 
the flux from positive relative velocities is 1.64 Jy km s~^, 
around twice as much as the flux from negative relative ve- 
locities of 0.91 Jy km s~^. The ionised gas does not show 
the same bias as the molecular gas, with fairly symmetric 
emission in both space and velocity. In addition, the extent 
of the ionised gas is much greater than that observ ed in the 
molecular gas (see Fig. 14)1 . As shown in Wiklind fc Henkel 
(|200ll ). however, the dust favours the north, in agreement 
with the molecular gas bias. 

Integrating the spectrum over the velocity range with 
observed emission, i.e. from -80 to 80 km s~^, we obtain a 
total CO flux of 2.77 Jy km s^^. The total flux is used to 
compute the total molecular hydrogen mass using the for- 
mula M(H2) = (1.22 X 1O''M0)D^ X ^co, where D is the 
distance measured in Mpc and ^co is the total CO(l-O) flux. 
This formula comes from using the standard CO to H2 con- 
version ratio iV(H2)//(CO) = 3 x 10^'' cm-^ where iV(H2) 
is the column density of H2 and /(CO) is the CO(l-O) inten- 
sity in K km s~^. We note that the actual conversion ratio 
for NGC 4550 is unknow n and can be expe cted to vary by a 
factor of at least two (e.g. lLerov et al.l l2007'). However, using 
the specified conversion ratio, the total H2 mass detected by 
our interferometric observations is 8.1 x 10^ ^Q- 

Using the same value of A^(H2)//(CO) and the same 
distance to NGC 4550, the single-dish measurement of 
IWikhnd fc Henkell l|200ll) gives a molecular mass of 1.4 x 10 
Mq, slightly less than twice the value we obtain. However, 
they integrate over a much larger velocity range (225 to 535 
km s~^), likely including much noise as positive emission 
(our velocity range is only 355 to 515 km s^^). Over a sim- 
ilar velocity range, our estimates for the molecular mass in 
NGC 4550 would be in suitable agreement, considering cal- 
ibration uncertainties (~ 20%). Reflecting the calibration 
uncertainty and the much larger uncertainty in the conver- 
sion ratio, we report the molecular mass in NGC 4550 as 
1 X 10^ 1^01 which is the value we will use for the rest of 
the paper. 



3 NUMERICAL SIMULATION 

While both the gas accretion and coplanar major merger 
scenarios seem to be able to produce at least the general 
properties of the two stellar discs in NGC 4550, we would 
also like to know how the gas in this galaxy came to ro- 
tate like the thicker disc. Gas rotating like the thick disc 
is actually in disagreement with the predictions of the gas 
accretion scenario, in which the thin disc is formed from 
the acquired cold gas (|Kennev et al.ll2000l ). Any remammg 
gas would then be expected to rotate like the thin disc. To 
see what the merger scenario predicts for the gas, we have 
performed a numerical simulation, described below. 
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Figure 3. Spectrum of the spatial region with observed CO(l-O) 
emission in NGC 4550. 
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Figure 4. CO(l-O) conto urs over the H/3 emission map of 
NGC 4550 from SAURON l lSarzi et al.ll200d) . The H/3 emission 
is clearly more extended than the CO. CO(l-O) contours are at 
0.25, 0.75, 1.25 and 1.75 Jy beam"! km s"!. 



3.1 Method and initial conditions 

Our merger simulation for NGC 4550 has been run with the 
same TREE-SPH technique as described in Di Matteo et al. 
(2007). The two merging galaxies are initially of the same 
type (Sb) and mass (each has a total mass 2.4 x 10^^ M© 
including dark matter). The bulges and halos are modelled 
as Plummer spheres, of masses 1.15 and 17.25x10^° Mq, and 
characteristic radii 1 and 12 kpc, respectively. The stellar 
discs have masses of 4.6 x 10^° M© and radial scale lengths 
of 5 kpc. The gas to stellar disc mass ratio is 0.2. 

The initial separation between the two galaxies is 100 
kpc. We choose an orbit for the galaxies such that they would 
have a parabolic encounter if they were both point masses. 
If the galaxies followed the parabolic orbit, the first peri- 
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Figure 5. Particle plots (x-y projection) of six snapshots of the simulation: the face-on views of the gas (left, blue) and stars (right, 
red) are plotted for each epoch, T=250 to 1500 Myr in steps of 250 Myr. The relative orbital motion is counter-clockwise in this picture, 
where the prograde galaxy is labeled "P" and the retrograde one "R". The merger is complete at the 3rd snapshot (T= 750 Myr), with 
some gaseous tidal dwarfs remaining thereafter. All boxes are 80 kpc wide. 



centre would be 8 kpc and their relative velocity at this 
point would be 707 km s~^. Of course, the model galaxies 
do not follow this keplerian orbit. Instead, given the strong 
dynamical friction, merging occurs rapidly. The orbital an- 
gular momentum is initially oriented in the positive z-axis 
Eis is the spin of one of the two galaxies (hereafter called the 
prograde galaxy). The other galaxy has negative spin (and 
will be called the retrograde galaxy). 

The total number of particles is 240 000, with 120 000 
in each galaxy, divided between 40 000 in the gas, 40 000 
in the stars and 40 000 in the dark matter halo. We use 
the Schmidt law and hybrid particles to take star formation 
into account as in Di Matteo et al. (2007). In the following, 
when we refer to the "gas component" , this means the sum 
of the gas particles and the very young stars formed since 
the merger. 



3.2 Analysis and Results 

Figure [S] shows the gas and stellar particles for six snapshots 
from 250 to 1500 Myr. The simulation was run until T=4000 
Myr (with a timestep of 0.5 Myr), but the evolution after 
1500 Myr is only secular relaxation of the merger remnant. 

With the spin of one galaxy aligned with the orbital mo- 
tion and the spin of the other anti-aligned, the interaction 
is strongly asymmetric despite the equal mass of the galax- 



ies. The prograde galaxy develops long tidal tails, while the 
retrograde is less heavily influenced by the interaction. As 
the merger evolves, the central regions become dominated 
by the stars from the retrograde galaxy as it remains more 
compact. However, the gas behaves differently due to its col- 
lisional nature. While two counter-rotating systems of stars 
can coexist, one of the gas systems must dominate. Initially, 
the gas from the prograde galaxy expands in tidal tails and 
even tidal dwarfs. However, this gas then falls back and set- 
tles in the disc of the merger remnant due to its dissipa- 
tive nature. Since the orbital angular momentum is posi- 
tive, the prograde gas increases its angular momentum and 
the retrograde gas loses some of its angular momentum in 
the encounter. This results in the gas rotating with posi- 
tive angular momentum after the gas from the two discs has 
interacted and fully settled down. 

Figure [6] displays both the edge-on gas density and ve- 
locity map for two time steps: it explicitly reveals how the 
gas settles into a thin disc while re-ordering its kinematics to 
prograde rotation. More quantitatively, Figure [7] shows the 
global rotation curves along an edge-on projection, splitting 
the rotational profiles according to the various components, 
the gas and stars in the two galaxies. The rotation curves 
make it very clear that the stars in the retrograde galaxy 
retain their direction of rotation, dominating the total ro- 
tation curve in the centre of the galaxy. The stars from the 
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Figure 6. The edge-on (x-z) stellar density (left), gas density contours (middle) and gas velocity field (right) are plotted for two times 
after the merger, T=1000 and 1250 Myr. It is easy to see that the settling of the gas occurs during this epoch: the gas density contours 
flatten and the gas velocity field becomes regular. The boxes are all 10 kpc wide, the wedges at the right of the velocity maps indicate 
the amplitude of the projected velocities in unit of 100 km s~^. All later snapshots are similar to the T=1250 Myr one, with a regular 
rotation velocity map for the gas. 



prograde galaxy show a lower amplitude of rotation, as they 
have been strongly heated in the encounter. The merger has 
thus created a disc galaxy with one dynamically hotter disc 
of stars, one dynamically cooler counter-rotating disc of stars 
and a gaseous component rotating like the thick disc. As the 
alignment of the orbital angular momentum with the pro- 
grade galaxy causes both the heating of this galaxy and the 
settling of the gas into prograde rotation, we predict that 
in such systems the gas will always be in corotation with 
the most perturbed stellar system, i.e. with the hotter disc. 
NGC 4550 clearly exhibits this phenomenon with both the 
molecular and ionised gas rotating like the thicker stellar 
disc. 

However, in the present simulation with perfectly copla- 
nar galax;ies, we note that there is not much heating in the 
perpendicular direction. The thickness of the stars with pos- 
itive angular momentum (those that would observationally 
be considered in the prograde disc) is comparable to the 
thickness of the stars with negative angular momentum. 
While the lack of a bona fide thicker disc limits the direct 
comparison to NGC 4550, we expect more resonant heating 
in the perpendicular direction if the interaction were not 
perfectly coplanar. Earlier simulations show that the an- 
gle between the discs' orientation and their orbit must be 
significa nt before the remnant takes on an elliptical mor- 
phology (jBournaud. Jog fc Combe3l2005l ). leaving room for 



disc thickening before disc destruction. We will explore the 
possible parameters of low-inclination mergers in a future 
work. 

As angular momentum exchange strongly drives the 
evolution of the merging galaxies, it is interesting to follow 
the angular momentum exchange between all components 
present, including the dark matter. Figure |S] shows clearly 
that the largest angular momentum exchanges occur during 
the violent merging epoch between 400 and 800 Myr. After- 
wards, all momenta evolve more slowly. The main feature is 
the loss of all the initially positive orbital angular momen- 
tum in each component, which is then transferred almost 
exclusively to the dark matter, except for a tiny fraction 
transferred to the internal spin of the gas and stars. This 
exchange is seen in the total global momentum (black full 
curve) in each panel of Figure [S] It is particularly intriguing 
to note that both dark haloes acquire a large prograde rota- 
tion, an effect even more marked for the halo of the prograde 
galaxy. 

Our simulation agrees with the results of Puerari & 
Pfenniger (2001) who first showed that a coplanar major 
merger between two disc galaxies could limit the heating ef- 
fect and create a NGC 4550 look-alike. In their simulation, 
they consider either a parabolic or a circular orbit, assum- 
ing that the relative energy of the two galaxies coming from 
infinity has already been absorbed by some outer matter. 
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Figure 7. Velocity evolution during the simulation. For the same epochs as Figure 6, the left panel shows the stellar density (solid black 
line), mean line-of-sight velocity (black dashed line) and velocity dispersion (dot-dashed red line) profiles. The central panel shows the 
mean stellar velocity profile of the resulting system again (solid black line), as well as those of the prograde (solid red line) and retrograde 
(solid blue line) galaxies individually. The right panel is as the central panel but for the gas velocities. All the profiles were extracted 
from a 4.5 kpc thick horizonthal slice taken from the edge-on projection shown in Figure 6. All boxes are 30 kpc wide. 



Both orbits avoid excessive heating, but only the parabolic 
orbit produces the spectacular counter-rotation observed in 
NGC 4550. However, their simulations are limited in how 
they investigate the gas- they include gas in either the pro- 
grade galaxy or the retrograde galaxy, but do not consider 
the case where gas is present in both galaxies as we have 
done. Of course, the final distribution of the gas within the 
merger remnant will depend on its initial distribution in the 
two progenitor galaxies. 

iDi Matteo et al.] l|2008l ) have also studied dynamical 
mechanisms to produce counter-rotating systems, through 
the merging of a non-rotating early-type galaxy with a gas- 
rich spiral. In these types of mergers, the orbital momen- 
tum is transferred to the rotation of the early-type galaxy 
and to the outer parts of the spiral, while the centre of the 
spiral (gas and stars) keeps its retrograde momentum. The 
counter rotation is then spatially marked (centre versus the 
outer parts), which is a different pattern than that observed 
in NGC 4550, where gas and stars are counter-rotating at 
the same location. The nature, amount and extent of the 
observed counter-rotation can thus help to determine the 
merging history of the system. 



4 STAR FORMATION 

The CO detected in NGC 4550 indicates the potential for 
star formation, as stars originate in cold molecular gas. 
Using the Toomre density criterion, Q, with a conserva- 
tively cold velocity dispersion of 6 km s^^ and an epicyclic 
frequ ency measur ed from the rotation of the ionised gas 
(Sar zi et al.ll200(^ . we calculate a threshold density (Q=l) 
of 285 Mq pc~'^. Determining a gas surface density for 
NGC 4550 from the CO data is difficult given the lack of 
resolution, but assuming a uniform-density, symmetric 5"- 
radius disc, we obtain a surface density of only 18.3 Mq pc~^ 
for the molecular gas, far below the threshold. Even decreas- 
ing the radius to 2", about the minimum believable extent 
of the molecular gas, gives a surface density of only 115 
Mq pc~^, still below the star formation threshold density. 
However, if the molecular gas is significantly more clumpy 
than a uniform density disc, than there may be isolated re- 
gions where the critical density is reached. We therefore do 
not immediately rule out the possibility of star formation in 
NGC 4550, instead looking for signs in other star formation 
tracers. 

Far-infrared (FIR) emission is often used as a tracer 
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Figure 8. Angular momentum evolution during the simulation. The top left panel shows the sum of all components, while the other 
three are split into the dark matter halo, the stars and the gas, as indicated above each panel. The wiggles in the total angular momentum 
line in the upper left box reflect noise in the simulation of order 2-3%. The colour scheme is as for middle and right panels of Figure 
7, i.e. red for the prograde galaxy and blue for the retrograde galaxy. The solid lines indicate the orbital angular momenta, while the 
dashed lines indicate the internal spin momenta. The angular momentum is in units of 2.3 X 10^^ Mq kpc km s^^; note that the scale 
is different in each plot. 



of star formation, as young OB stars efficiently heat their 
surrounding dust. In late-type galaxies where the ultravio- 
let (UV) and visible radiation is dominated by young stars, 
the FIR emission is a very good tracer of star formation. 
Its use in earlier-type galaxies, however, is more disputable, 
as the radiation field of old hot stars should also contribute 
to the dust heating. Yet the relatio ns between other tr acers 



of star formation such as Ha fe.g. iKewley et al.|[20o3 ) and 



Ley ' 

radio continuum (e.g. iGavazzi et al.lll986l ) hold puzzlingly 
constant with morphology (at least within spirals), indicat- 
ing that it is dust proximate to star forming regions which 
dominates the FIR emission. Converting the IRAS 60 and 
100 jxm fluxes i nto a star formation rate (SFR) following the 
prescription of iKewlev et al.l (|2002| ). we derive a FIR-based 
SFR of K. 0.017 Mq yr"^ for NGC 4550. 

As many early-type galaxies without star formation also 
emit in the FIR (an active galactic nucleus or infrared cirrus 
can also contribute) , we do not want to rely on the FIR emis- 
sion alone to ascertain ongoing star formation in NGC 4550. 
Instead, we look to the UV, where young stars are particu- 
larly bright. The first diagnostic we consider is the NUV-V 
colour. Using optical data from the MDM 1.3m telescope 
(Falcon-Barroso et al. in preparation) and GALEX UV data 




20 30 40 

Distance (arcsec) 

Figure 9. NUV and V surface brightness profiles and NUV-V 
colour along the major axis in NGC 4550. The vertical dashed 
line indicates the effective radius while the horizontal dashed line 
indicates the blue limit to the colours the uv-upturn phenomenon 
can produce, NUV-V=5.2. 



( Jeong et al. in preparation) , we have produced a plot of the 
NUV-V colour along the major axis of NGC 4550 (Fig. [Q]). 
NGC 4550 has a central NUV-V colour of around 5, becom- 



Molecular gas and star formation in the red-sequence counter-rotating disc galaxy NGC 4550 9 



ing bluer with increasing radius. While old stars can also 
emit in the UV, pro ducing the UV-uptu rn phenomenon seen 
in old galaxies (e.g. lBurstein et al.lll988l . see O'Connell 1999 
for a review), t he bluest NUV -V colour this eflect produces 
is NUV-V=5.2 (|Yi et al.ll2005h . Thus, the NUV-V colour in 
NGC 4550 is more likely from young stars. Furthermore, the 
decreasing gradient (bluer colour with increasing radius) is 
not seen in UV-upturn galaxies, reinforcing the idea that the 
NUV emission comes from young stars. 

To inves tigate this possibili ty, we have used the method 
described in Ijeong et al.l l|2007l ) to fit a two-stage star for- 
mation history to the UV-optical colours of NGC 4550. An 
old population is fixed at an age of 10 Gyr with a composite 
metallicity, while a younger component at solar metallic- 
ity is allowed to vary in age (0.001 Gyr < Iyc < 10 Gyr) 
and mass fraction (10"'* < /yc < 1). Figure [TOl shows the 
results of these fits, pixel by pixel. The central region ap- 
pears to host a young stellar population (fitted age around 
100 Myr), although the mass fraction is very low at only 
0.01%. To increase the signal-to-noise and thus the reliabil- 
ity of our model fits, we computed integrated colours for the 
central region where the CO emission is detected, using an 
ellipse with a major axis of 5", a minor axis of 3" and a 
position angle of 0°. The colours from this aperture give a 
best-fit age of 280 Myr with 5.89 x lO'^ Mq of young stars. 
The fit to the spectral energy distribution (SED) of this cen- 
tral region of NGC 4550 is shown in Figure 1111 As seen in 
the figure, NGC 4550 is slightly bluer in NUV-V than the 
strong UV-upturn galaxy NGC 4552 and much bluer than 
the intermediate age compact elliptical M32, indicating that 
young stars must be responsible for the blue NUV-V colour 
of NGC 4550. Figure [TT] also displays an inset of the chi- 
square contours for the fit. As can be seen, an age-mass de- 
generacy is present, with larger mass of slightly older stars 
(but only up to 500 Myr) or a smaller mass of even younger 
stars also able to fit the colours. 

Absorption linestrengths also give an indication of 
a young or intermediate age population in NGC 4550. 
While absorption linestrengths are not sensitive to the very 
youngest stars, they can easily reveal stellar populations 
1 Gyr or older. The age-sensitive linestrength H/3 (measured 
on the Lick/IDS system) is relatively high in NGC 4550, at 
a value of 1.99 A integrated over an effective radius (Re ) 
and 2.14 A integrated over Re/8 l|Kuntschner et al.l [20061 ). 
reaching a value of 2.2 A in the very centre (Maier et al., in 
preparation). Observations with the Multi-Pupil Fiber Spec- 
trograph (MPFS) measured lower H/3 linestrengths, how- 
ever these were not c orrected for the contaminati on from 
the H/3 emission line (|Afanasiev fc Sirchenkoll2002l ). These 
linestrength values are hard to reproduce with only old (> 5 
Gyr) stars and more likely indicate a young or intermedi- 
ate age population, perhaps dating from the time of the 
merger. As with the UV, there is again an age-mass fraction 
degeneracy, with a smaller fraction of young stars producing 
an analogous effect to a larger fraction of intermediate-age 
stars. 

Unfortunately, neither the UV nor the absorption 
linestrengths can conclude anything about ongoing star for- 
mation in NGC 4550, while the FIR only gives an upper 
limit. To investigate ongoing star formation further, we look 
at optical emission line ratios and polyaromatic hydrocarbon 
(PAH) emission. Particularly low ratios of [OIII] to H/3 emis- 



sion lines [log([0 IIl]/H/3) < —0.2] indicate curren t star for- 
matio n (< 20 Myr) as the source of ionisation ||Ho et al.l 
Il997h . Data from SAURON indicate th at log([0 III] /H/3 ) 
is between 0.15 and 0.5 for NGC 4550 l|Sarzi et al.ll2006l ). 
While gas could be star forming at these high values, it is 
more probable that shocks, post-asymptotic giant branch 
stars or an AGN are the dominant source of ionisation 
(Sarzi et al. in preparation). The greater physical extent of 
the ionised gas compared to the molecular gas (see Fig. |3| 
also supports the idea that star formation is not the pri- 
mary ionisation source. Similarly , although the PAH emis- 
sion detected in NGC 4550 ( Brcs san et al.ll2006D IS sugges- 
tive of star formation, the PAH spectrum lacks the dis- 
tinctive line ratios seen in star-forming regions. This dif- 
ferent spe ctrum may just refle ct a very low star forma- 
tion rate (|Galliano et al.l |2008| ). but it may also indicate 
anoth er source of PAH excitation altogether (Smith et al.l 
l2007f) . Thus, with the available data, we cannot say more 
than that the star formation rate in NGC 4550 is less than 
about 0.02 Mq yr~^, as indicated by the FIR emission. 

At this rate, the amount of molecular gas we detect 
in NGC 4550 can fuel star formation for up to 350 Myr. 
The molecular gas is likely related to the young population 
(280 Myr) found in the UV, although whether the star for- 
mation has proceeded continuously or in a more episodic 
fashion since then is impossible to say. Depending on the in- 
terpretation of the linestrengths, the star formation period 
may be of an even longer duration (extending back to more 
than 1 Gyr) or the moderately high linestrengths may reflect 
a stellar population originating in the merger that formed 
NGC 4550. 



5 CONCLUSIONS 

We detect a very small amount (1 x lO'^ ^q) of molecu- 
lar gas in the centre of the counter-rotating red disc galaxy 
NGC 4550. The CO(l-O) emission is limited to the central 
750 pc (10") and is asymmetrically distributed, stronger at 
positive relative velocities north of the galaxy centre. The 
molecular gas co-rotates with the thick stellar disc and the 
ionised gas, counter-rotating with respect to the thin stellar 
disc. 

Our simulation of the merger of two counter-rotating 
coplanar disc galaxies shows that the main features of 
NGC 4550 can naturally be explained with such a scenario. 
The interaction heats the prograde disc more than its ret- 
rograde companion and the gas component ends up aligned 
with the total angular momentum (dominated by the orbital 
angular momentum), and thus with the prograde disc. Thus, 
the gas is predicted to always rotate like the dynamically 
hotter disc, as observed in NGC 4550. As the gas accretion 
scenario does not provide a natural explanation for the coro- 
tation of the gas and thick disc, the merger scenario appears 
the more likely. 

Both the UV and optical-linestrength data indicate that 
NGC 4550 cannot be made up of a purely old stellar popu- 
lation. The best-fit two-population model to the UV-optical 
photometry in the region with observed molecular gas gives 
a young population of 280 Myr and mass of 5.89 X 10^ Mq. 
The optical H/3 linestrength value requires a population at 
least younger than 5 Gyr. Ongoing star formation in the 
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Figure 10. Left: Two-component fit map of the age of the young component (in Gyr). Right: Two-component fit map of the mass 
fraction of the young component. Note that the mass fraction is relevant to each individual pixel, thus a certain mass fraction in the 
central pixels represents a larger mass of young stars than the same mass fraction in outer pixels. 
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Figure 11. Integrated V, NUV and FUV fluxes for the central region of NGC 4550. The central region is defined by an ellipse of 5" x3", 
PA=0°. The spectral energy distributions (SEDs) of the best-fit young and old components are shown, along with the SED of the sum 
of these components. We have also shown the integrated colours of NGC 4552 (strong UV-upturn galaxy) and M32 (intermediate-aged 

compact elliptical) normalized in the V-band for comparison. Inset: contours of the two-component stellar population fit to the central 
integrated colours of NGC 4550, showing the age-mass fraction degeneracy. The best fit is marked with an 'x'. 



observed molecular gas is possible if there axe locally dense 
regions that exceed the critical density. The FIR emission 
gives an upper limit on the current star formation rate of 
0.02 Mq yr^^. This low star formation rate combined with 
the small amount of molecular gas present suggest that we 
are either witnessing a weak period in an extended bursty 
star formation episode, or a more continuous very low-level 
star formation episode that will last about another 350 Myr. 
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